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It is well known that conventional pairing fluctuation theory at the Hartree level leads to a normal state pseu- 
dogap in the fermionic spectrum. This pseudogap arises from bosonic degrees of freedom which appear at 
T* , slightly above the superconducting transition temperature T c . Our goal is to extend this Hartree approxi- 
mated scheme to arrive at a generalized mean field theory of pseudogapped superconductors for all temperatures 
T. While an equivalent approach to the pseudogap has been derived elsewhere using a more formal Green's 
function decoupling scheme, in this paper we re-interpret this mean field theory and BCS theory as well, and 
demonstrate how they naturally relate to ideal Bose gas condensation. Here we recast the Hartree approximated 
Ginzburg-Landau self consistent equations in a T-matrix form, from which we infer that the condition that the 
pair propagator have zero chemical potential at all T < T c is equivalent to a statement that the fermionic exci- 
tation gap satisfies the usual BCS gap equation. This recasting makes it possible to consider arbitrarily strong 
attractive coupling, where bosonic degrees of freedom appear at T* considerably above T c . The implications 
for transport both above and below T c are discussed. Below T c we find two types of contributions. Those as- 
sociated with fermionic excitations have the usual BCS functional form. That they depend on the magnitude of 
the excitation gap, nevertheless, leads to rather atypical transport properties in the strong coupling limit, where 
this gap (as distinct from the order parameter) is virtually T-independent. In addition, there are bosonic terms 
arising from non-condensed pairs whose transport properties are shown here to be reasonably well described by 
an effective time-dependent Ginzburg-Landau theory. 



I. INTRODUCTION 

In this paper we address the physics of the pseudogap state 
on the basis of a beyond-BCS mean field theory. Our approach 
represents a natural extension to arbitrarily strong coupling 
constants, g, of time-dependent Ginzburg-Landau (TDGL) 
theory when the quartic terms are treated at the Hartree ap- 
proximation level. We also explore the broken symmetry 
phase, incorporating the effects of a Hartree approximation 
to pairing (or equivalently, excitation-gap) fluctuations for 
T < T c in a fashion which is necessarily compatible with their 
description above T c . A central assumption in our work is that 
the ground state belongs to the generalized BCS-like theory 
introduced by Leggett 1 , and further analyzed by Nozieres and 
Schmitt-Rink 2 , and others. In this way, one arrives at a simple 
re-interpretation of BCS theory as a limiting case of a more 
general theory of ideal Bose gas condensation. In contrast to 
true Bose systems, here, however, fermionic degrees of free- 
dom are never absent: the condition fi pa ir = 0, for T < T c is 
equivalent to the constraint imposed by the BCS gap equation 
on the excitation gap A(T). 

Several key points need to be made. The observation of 
a pseudogap (T > T c ) deriving from Hartree approximated 
treatment of "pairing" fluctuations dates back to the early 70's, 
with the work of Abeles 3 , PattorA 5 and Hassing and Wilkins 6 . 
Pairing fluctuations can be alternatively viewed as fluctuations 
in the magnitude of the fermionic excitation gap, that is "gap 
fluctuations". Most importantly, at the Hartree level, these are 
to be distinguished from order parameter fluctuations, as will 
be discussed in more detail below. Indeed, bosonic degrees 
of freedom, associated with fermion pairs, and the existence 
of a fermionic excitation gap are two sides of the same coin. 
The gap or pseudogap in the fermionic excitation spectrum 



reflects the fact that extra energy is needed to break bosons, 
or fermion pairs, apart. When approached from above, T c is 
suppressed (relative to its strict BCS counterpart) by the exis- 
tence of this normal state pseudogap, as a consequence of the 
lowering of the fermionic density of states near Ep. When 
approached from below, one readily sees that the superfluid 
density will not vanish at T* , where the fermionic gap disap- 
pears, but instead at a lower temperature T c , as a consequence 
of additional (beyond BCS) bosonic excitations of the con- 
densate. 

In conventional {i.e., long coherence length, £) three dimen- 
sional superconductors, these non-condensed pairs are virtu- 
ally undetectable. However, there are two mechanisms for en- 
hancing their effects: reduced dimensionality, often with dis- 
order [see the recent review 7 by Varlamov and Larkin which 
addresses Gaussian pair fluctuations] and increased attractive 
coupling constant 2 *^ associated with short £. If by "the 
fluctuation regime" we mean the temperature range where 
bosonic degrees of freedom are present, then in the second 
mechanism one is effectively extending this regime up to T*, 
which could be as much as an order of magnitude higher 
than T c . In the vicinity of Bose condensation there are di- 
vergences in various transport properties, but even away from 
T c , closer to T* , bosonic contributions may dominate over 
their fermionic counterparts, leading to highly unconventional 
transport characteristics, as will be outlined here. 

The Hartree approximation&I2iii and its T-matrix 
extension 4 ' 5 is a previously well established route for ob- 
taining a pseudogap state in homogeneous superconductors 3 . 
Throughout this paper we call the latter "extended Hartree 
theory". Nevertheless, there is considerable emphasis in 
the current literature on approaches to the pseudogap based 
on an alternative phase fluctuation scenario 12 . The contrast 
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between these two schemes should be emphasized. In the 
former case one is dealing with non-condensed bosons, 
which are associated with fluctuations in the fermionic 
excitation gap A(r, t), (say, measured relative to the spatially 
uniform or mean field order parameter). In the latter case, 
one is dealing with fluctuations in the phase of the complex 
order parameter A sc (r,t). At a beyond-BCS level these 
two channels are distinguishable (A ^ A sc ), but, as in 
any mean field theory, the order parameter is taken to be 
spatially uniform. By contrast, the phase fluctuation scheme 
is based on order parameter fluctuations around the strict 
BCS state (A = A sc ). Generally, order parameter fluctuation 
approaches are appropriate when the separation between the 
mean field onset temperature T* and the actual instability 
temperature T c is small. If this temperature difference is not 
small, it is more appropriate to derive an improved mean field 
theory first 13 and then append fluctuation effects, if they are 
called for. This is the philosophy of the present approach in 
which at a mean field level, the separation between T* and T c 
can be arbitrarily tuned via the size of the attractive coupling 
constant. 

Finally, we emphasize a contrast between the present ap- 
proach and previous work by Nozieres and Schmitt-Rink^and 
by Randeria and co-workers^. Both groups studied the ef- 
fects of pre-formed pairs, or non-condensed bosons, presum- 
ing a BCS-like ground state with arbitrary attractive coupling 
g. We will assume this ground state as well. These other 
studies were conducted at the level of Gaussian pair fluctua- 
tions, so that pseudogap effects were absent in the pair prop- 
agator. Moreover, an approximate form for the number equa- 
tion was assumed as well 14 . Essential to the physics we will be 
discussing, is a pseudogap in the fermionic spectrum, which 
helps to decouple the bosonic and fermionic degrees of free- 
dom. This pseudogap leads to relatively long-lived bosons, 
without the need to invoke unphysically strong coupling. A 
gap in the fermionic spectrum inhibits the decay of bosons 
into fermions, simply because there are very few low energy 
fermionic states to decay into. As a consequence, bosonic de- 
grees of freedom are evident over a much wider temperature 
range, above and below T c , once g is beyond the weak cou- 
pling regime. 

In this paper, we make no contact with experiments in high 
temperature superconductors. The goal here is to demonstrate 
the simple physics, as well as the genesis of pair fluctuation 
approaches, as distinguished from phase fluctuation schemes. 
It is useful to emphasize that we are focusing in the present 
paper on the nature of the superconductivity (below T c ) and 
its implications above T c . While we presume a generalized 
BCS ground state (albeit, with arbitrary g and self consistent 
fermionic chemical potential), other ground states have been 
contemplated in the cuprate literature which accommodate the 
physics of the Mott insulator to varying degrees. The strongest 
support for the relevance of our viewpoint (which effectively 
sidesteps Mott effects) is the anomalously short coherence 
length £. This suggests a breakdown of strict BCS physics, 
which, at the very least, needs to be understood and character- 
ized, on its own, independent of Mott effects. To support this 
breakdown is the highly non-BCS temperature dependence of 



the measuredi£ii& excitation gap A - not so different from 
that shown in Fig. [3] below. Although transport calculations 
and interpretations of transport data are often based on a BCS 
ground state, rather little attention has been paid in the past 
to the anomalous behavior of A(T) and its implications for 
transport. In this paper we address this omission. 



II. ABOVE T c AT WEAK COUPLING 

A. Overview of Hartree approximated Ginzburg-Landau 
theory 

The Ginzburg-Landau (GL) free energy functional in mo- 
mentum space is given by£ 

P Q 

+ iE 6 *eA%*Q 2 %%+Q J -«3 . (i) 
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where are the Fourier components of the order parameter 
*(r, t), Q = (i£l n , q), e = a = £i is the GL co- 

herence length, T* is the critical temperature when feedback 
effects from the quartic term are neglected, j3 — l/T (ks is 
set to 1) and N(0) is the density of states at the Fermi level 
in the normal state. We approximate the quartic term so that 
only paired terms are included in the last addend of Eq. Q 
leading to 
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That there is no factor of 2 in the first term on the right- 
hand side of the above expression reflects the fact that we use 
Hartree rather than the Hartree-Fock approximation. As found 
elsewhere 6 , fo y = bQ i q j can be approximated by fcofoiofojO 
where b = [A^(0)V/7T 2 ]|C(3)- To further simplify the quar- 
tic term, we apply the mean field approximation, writing 
* q o| 2 = (|* q o| 2 ) + (5|* q0 | 2 and neglecting in Eq. (0 terms 
of order (<5|>I< q0 | 2 ) 2 . This leads to 
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The contribution ( | 1 2 ) is determined self-consistently via 
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when we replace the quartic term in Eq. Q by Eq. Q. It 
follows that 
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If we sumEq. @ over q and identify X) q (l*qo| 2 } with /? 2 A 2 , 
we obtain a self-consistency equation for the energy "gap" (or 
pseudogap) A above T c 
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Note that the critical temperature is renormalized downward 
with respect to T* and satisfies 
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B. Introduction to T-matrix : T ^ T c , small A(T C ) 
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Equations |7}, (|8}, and l|9} are the central equations de- 
rived from the Hartree approximated Ginzburg-Landau (GL) 
scheme. They describe how the excitation gap A(T) and the 
quantity /2 pa , r behaves above, but near T c . We now rewrite 
these equations using a T-matrix approach. 

A central quantity in T-matrix schemes is the pair suscepti- 
bility. Here we take this function to be of the form^ 



x(Q) = J2 G °(Q ~ K )°( K )vl 
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where Go is the bare Green's function and G the full or 
dressed Green's function which depends on the self energy 
T,{K ) given by 
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Here and throughout we use the convention J^q = 
T J2in S q - l n tne a b° ve expressions tp^ represents a gen- 
eralized (for example) s- or d-wave function symmetry factor. 
While there are two other T-matrix approaches in the litera- 
ture, one2£ in which \ is related to GqGq, and one which 
is even more widely usedii in which \ is related to GG, 
only the GGo scheme is simply connected to the Hartree-GL 
approach^ 

To compare with GL theory we expand these equations to 
first order in the self energy correction 18 . The T-matrix can be 
written in terms of the attractive coupling constant g as 
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Defining 



A 2 = -£i(Q) 
Q 



we arrive at (see AppendixlATl 

E(A) « -G (-A)AV 2 
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The results of Appendix IA1 can be used to derive a self con- 
sistency condition on A 2 in terms of the quantity Sx(Q) (first 
order in £), which satisfies 



implying that 



6 X (0) 



Sx(0) = -M/3A) 2 , 



b f d 3 q 



(16) 



N(0)Tj (2tt)3 e + gq2-5 X (0)/N(0) 



(17) 

which coincides with the condition derived earlier in Eq. JSJ. 
Finally, defining 

»pair(l ) - N{Q)m ~ £+ A(0) ( ^ 

we may interpret the vanishing of p, pa i r as the condition that 

at T = T r 
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It follows from Eqs. {T^, {T^J and © that 
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so the above condition for T c is in agreement with that found 
earlier [Eq. (|9j] and the effect of a finite A(T C ) (self-energy 
correction) is a shift downward in the critical temperature rel- 
ative to its value (given by T*) in the A(T C ) = limit. 



III. BELOW T c : WEAK COUPLING 

A. Hartree approximated GL theory for small A(T C ), T ^ T c 

The left hand side of Eq. may be interpreted as the num- 
ber density of bosons of momentum q. Since p,p a ir{T c ) = 0, 
the q — level becomes macroscopically occupied once the 
system enters the superconducting region, at T = T c . To sup- 
port this assertion we investigate the behavior of p, pa i r (T) for 
T « T~. We separate out the q = term in Eq. (J6J and write 
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Eq. (|8j leads to another constraint on A, which yields 
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where c = bo/[N(0)V}. Differentiating both of the above 
equations with respect to T one obtains the following expres- 
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The numerator of Eq. d23l > is negative in the vicinity of T c , 
A ' a i t 1 me fi]- s t term in the denominator di- 



since 



< fc. At T c . 



verges ( v - is a finite number in the thermodynamic limit), 



and as also found elsewhere 



6 rfjJp 



dT 



(T = T c ). Since 



Jipair cannot be positive (that would make the right hand side 
of Eq. negative), and its derivative is negative or zero, we 
conclude that fL pa ir must remain zero in the vicinity of, but 
below T c , where the GL description is applicable: 



This implies, following Eq. l|8), 
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This result was previously obtained in Ref. El It is important 
because it shows that when the system reaches T c the excita- 
tion gap A(T) assumes the BCS or mean field value. 

We notice strong analogies with Bose-Einstein condensa- 
tion (BEC) in an ideal Bose gas. Here A 2 plays the role of 
N, the total number of bosons, which below T c contains two 
components, one associated with the condensate A 2 =0 and the 
other with the non-condensed pairs A 2 ^. The latter are like 
the excited states of the BEC system. We write 
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It follows from Eqs. 0, and <EU that 



/3 2 A 2 



d 3 q 



(2tt) 3 N VT Hq 2 



(26) 



(27) 



Just as in the ideal Bose gas problem, A 2 (T) is constrained 
(via the pair chemical potential condition), A g ^o (T) is con- 
strained through the self-consistent Hartree condition and thus 
one may deduce the condensate weight, or superconducting 
order parameter A 9=0 {T). 



B. Behavior near T = where A(T) is no longer small 

A useful observation can be made at this time. We have 
just shown that in the Hartree theory, A(T) assumes the BCS 
value in the vicinity of, but below T c . One expects on very 



general grounds that sufficiently far away from T c , pair "fluc- 
tuation" effects are irrelevant and that the system is described 
by strict BCS theory. (In this regime A(T) is no longer a 
small parameter). Thus we may infer that everywhere below 

T c 

A(T) - A BCS (T), T < T c , (28) 

so that the excitation gap is given by the BCS value. 10 What 
is different from strict BCS theory, however, is that 

A(T C ) + (29) 

This is the sole effect of pair fluctuations below T c . Never- 
theless it has important consequences, because it reflects the 
presence of non-condensed bosons below T c , which, in turn, 
mirror their normal state counterparts. 

C. T-matrix scheme below T c 

We now show that Eq. J2 81 is connected to the ideal Bose 
gas condition at all temperatures below T c . This follows from 
the analysis in AppendixlBlin which it is shown that the BCS 
gap equation is associated with a divergence of the T-matrix 
defined in Eq. dl 2i . at Q = for all temperatures below T c . 
Thus t~g(0) = implies 

Hpair{T) = Q, T<T C . (30) 

In order to satisfy this gap equation, together with an indepen- 
dent constraint on the number of finite momentum pairs one 
must incorporate a broken symmetry A sc ^ 0, which we now 
reformulate within our T-matrix theory. Below T c 
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can be decomposed into two contributions: the "pseudogap 
and superconducting components via 
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By defining 



(32) 

(33) 
(34) 

(35) 



and using the divergence of t(Q) at Q = to evaluate 
Eq. (13 II . we retrieve the BCS-like self-energy 

E(K) « -G (-K)AV 2 (36) 

with 

A 2 = A 2 C + A 2 fl (37) 

For small A(T C ) and T w T c , Eqs. «|30j,(|35j and 03 are man- 
ifestly equivalent to their GL counterparts (with fi pa i r corre- 
sponding to fipaw 19 , A sc to A q=0 and A pg to A q ^ ); they 
give a complete 9 description of the system by determining the 
full excitation gap A, order parameter A sc and the amplitude 
of propagating pairs (the pseudogap) A pg . 
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IV. EXTENDED HARTREE APPROXIMATION: 
ARBITRARILY STRONG COUPLING g 

We have demonstrated above that there is a simple T- 
matrix scheme involving the pair susceptibility (GGo) which 
is equivalent to Hartree approximated GL theory both above 
and below T c . (Moreover, there is a natural extension down 
to T = 0). This equivalence has been proven provided we 
restrict ourselves to small A(T C ), where GL approaches are 
applicable (see Fig.Q 

There is no reason, however, that the T-matrix scheme can- 
not be considered for arbitrary coupling constant g, or equiva- 
lently large A(T C ) where conventional GL theory is no longer 
appropriate. In this regime the separation between T* and 
T c can be extremely large. In this way, a pseudogap in the 
fermionic spectrum occurs at a temperature much higher than 
the Bose condensation temperature T c (see Fig. [2). As in con- 
ventional GL theory, this pseudogap reflects bosonic degrees 
of freedom. Once bosons are meta-stable it takes a finite exci- 
tation energy to create fermions from them. 

All the appropriate T-matrix equations have been presented 
above. [See Eqs. (go) and d32YQ7ll, Note that the only tech- 
nical difference between the cases of weak and strong cou- 
pling is in the details of the expression for the T-matrix it- 
self (see Appendix IBV There is, however, an important addi- 
tional constraint which needs to be appended, on the fermionic 
chemical potential p. 

The resulting equations greatly simplify at and below T c be- 
cause of the vanishing of the pair chemical potential. We sum- 
marize the above discussion by presenting expressions (appro- 
priate to the superconducting state) for the gap and chemical 
potential p 

ff -i + £izM^. =0 (38) 

k k 



n = 2^[f(E k )+i£(l-2f(E k ))] (39) 

k 

The former follows from the vanishing of t~g(0) in Eq. J32t 
and the latter from n = 2 J^k G(K). The quantity is the 
coherence factor — |(1 — ek/^k) with ek = k 2 / (2m) — p 
and is the fermionic excitation energy which depends on 
the magnitude of the superconducting order parameter A sc 
and the pseudogap energy scale A pg : 

E k = +A2(k) (40) 

where A 2 (k) = A 2 g (k) + A 2 c (k) = A 2 ^ 2 .. The important 
quantity A pg (T) here is deduced following Eq. (j35j. 

At a physical level the parameter which should be asso- 
ciated with strong or weak coupling is T* /T c . However, 
T* is frequently difficult to quantify. An alternative param- 
eter for characterizing the strength of the coupling is a = 
A(T C )/A(T = 0). This is more readily accessed experimen- 
tally. It should be seen from Fig. ^ that the BCS limit cor- 
responds to A(T c )/A(0) w 0, whereas (following Fig. [2} in 
the strong coupling regime A(T c )/A(0) » 1. Figure ^indi- 
cates how the temperature dependence of the gap varies with 
the strength of the coupling, or alternatively with a. This is 
an important plot because it epitomizes, perhaps, more than 
any other how dramatic are the differences from strict BCS 
theory. The challenge, then, is to accommodate the results 
in this figure into transport and other calculations within the 
superconducting phase. We thus turn to dynamical effects as- 
sociated with non-condensed bosons. 
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FIG. 3: Implications of the ideal Bose condensation condition 
(impair ~ 0, when T < T c ) for the excitation gaps. Each of the 
curves corresponds to different coupling constant strengths. 



tion, because pairing fluctuation bosons are essentially free; 
they interact with fermions but not directly with each other. 
Once there are boson-boson interactions (beyond mean field 
theory), then necessarily there is important coupling to order 
parameter fluctuations. 

Fermonic degrees of freedom also contribute to transport 
and thermodynamics, but there is no simple phenomenol- 
ogy (or counterpart of TDGL) for addressing these terms. 
When the boson contributions to transport are small (say, 
as in the thermal conductivity), the fermionic terms cannot 
be neglected, and these have to be computed diagrammati- 
cally, as will be discussed in more detail below. In other in- 
stances the bosonic contributions are singular 7 , or nearly so, 
in the vicinity of T c and TDGL calculations are appropriate 
for T C <T < T*, as well as for T < T c . In the remainder of 
this section we focus only on the bosonic terms. 



B. TDGL above and below T c 



V. EFFECTIVE TDGL THEORY FOR STRONG 
COUPLING LIMIT 

A. Overview 

In this section we investigate the applicability of an effec- 
tive time-dependent Ginzburg-Landau theory as a basis for 
characterizing the dynamics of our excitation-gap or pairing 
fluctuations. This discussion is relevant to transport calcula- 
tions associated with bosonic degrees of freedom. Quite gen- 
erally, we may view TDGL as a generic equation of motion 
for describing bosons interacting with a reservoir of fermionic 
pairs. This equation of motion can be associated with the long 
wavelength, low frequency behavior of the pair propagator or 
T-matrix, and it has a validity both above and below T c , as 
long as bosonic degrees of freedom are relatively stable. 

In the presence of a pseudogap, the situation is complicated 
by the fact that the order parameter is generally different from 
the excitation gap. This is illustrated clearly in Figs.^ an d 
Thus, below T c , one has to take care to distinguish be- 
tween the fluctuations in each of the two channels. By con- 
trast, above T c bosonic degrees of freedom are uniquely as- 
sociated with fluctuations in the fermionic excitation gap. To 
address order parameter fluctuations in the presence of a fi- 
nite gap A at T c has been done elsewhere^. It is conceptually 
straightforward, but difficult to implement. Such fluctuations 
also alter critical exponents within a generally narrow "critical 
fluctuation" regime. 

We will not consider these order parameter fluctuations 
here, but rather presume that A sc is spatially uniform, as in a 
mean field theoretic approach. The bosons of interest to us are 
non-condensed pairs which in turn lead to fluctuations in the 
fermionic excitation gap. A variation in the number of pairs 
leads to a variation in the excitation gap A, simply because 
one has to pull apart the constituents of the pairs to create 
fermions. Moreover, our mean field theory (associated with 
a BCS-like ground state 1 ) provides a considerable simplifica- 



At a microscopic level, one can expand the T-matrix for 
non-condensed pairs at small q, fi 



C(q,^) 



n 2 + (i + 



Q 

2M 



(41) 

In order to be consistent with the linear dynamics of TDGL, 
the quadratic terms in fl are neglected in what follows. The 
coefficients in Eq. J4U are in general T dependent, although, 
significantly below T*, the most important temperature de- 
pendence is associated with fi pa ir, which is zero at and below 
T c and negative above T c . The temperature T* corresponds 
to the onset of a resonant structure in the T-matrix. This onset 
temperature is reliably computed without including pseudo- 
gap effects- thus, at the Gaussian approximation level. Reso- 
nance effects -™ (reflecting the initial formation of meta-stable 
pairs) enter via the ratio a' Q /a,Q, which (at fixed T near T*) can 
be increased by increasing the coupling constant g. The larger 
is this parameter the more pronounced are pair resonances. 
Stated alternatively, as the coupling constant g increases the 
propagating component a' becomes increasingly more impor- 
tant than the diffusive term Oq. Indeed, similar observations 
were made by Randeria and co-workers 8 . 

We now proceed to the more detailed analysis of the coef- 
ficient Oq, slightly above and at all T below T c . Well above 
T c (but below T*) more detailed numerical calculations are 
requireoSS, and these demonstrate that Oq increases with 
T as pseudogap effects diminish. At our leading order ap- 
proximation as in Eq. ( I36> , we do not distinguish between 
lifetimes associated with the condensed and non-condensed 
bosons. For the purposes of computing T c and the various en- 
ergy gaps below T c , this has been shown to be a reasonable 
approximation^^, but it clearly misses important physics as- 
sociated with the onset of true phase coherence^!. In this ap- 
proximation the finite momentum pairs are extremely long- 
lived (see Eqs. JB 1 b and dCll l). Consequently, Oq is zero. For- 
mally, this results comes from the fact that the imaginary part 
of the pair susceptibility T(0, il) has a higher power than lin- 
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ear dependence on fl. If, instead we introduce (as in Appendix 
ICl a finite lifetime to distinguish the contributions from con- 
densed and non-condensed bosons 

m) = ^d + A >1 (42) 
lu + e k lj + £k + ij 

we find a non- vanishing 22 TDGL coefficient q'q above T c . 

< = » <43) 

Below T c , however, the existence of a condensate leads to 
<2q = 0, even for the more general self energy of Eq. i42l . 
These arguments are presented in Appendixlcl However, the 
presence of even a small amount of disorder is sufficient to 
restore a linear in frequency imaginary term in the TDGL ex- 
pansion of the T-matrix. These observations were made by 
Chen and Schrieffer— . In Fig.|4] we present a figure from their 
work which illustrates how the frequency dependent contribu- 
tions to the T-matrix evolve with impurity concentration, for 
different scattering strengths. 

In this way one establishes an effective TDGL description 
for the non-condensed boson dynamics both above and below 
T c . It should be stressed, however, that the character of this 
theory changes on either side of T c . 

The observation that the pairs live very much longer than 
anticipated by, say, Gaussian level calculations is a conse- 
quence of the fermionic pseudogap 22 - 23 . In this way, the diffu- 
sive component (i.e., the parameter o!q ) remains small for an 
extended range of temperatures above T c , and becomes even 
smaller (impurity limited) below T c . This reflects the fact that 
as the fermions acquire a larger gap, the bosons live longer and 
the two degrees of freedom become progressively more dis- 
tinct. These same effects are underlined by our earlier obser- 
vation that T c , as distinguished from T*, must be computed by 
including the feedback effects of the fermionic pseudogap on 
the bosonic propagator. This is precisely reflected in Eqs. (I38i 
and d39> which, together with Eq. fl37i must be solved to de- 
termine T c . Because of the possibility of a large separation 
between T* and T c , (as g is progressively increased) one may 
extend the simple dynamics associated with TDGL theory to 
describe bosonic transport at temperatures below T*, not just 
those limited to the immediate vicinity of Base condensation, 
T c . This provides some microscopic support for a recent phe- 
nomenological approach 26 which addresses Nernst and other 
normal state transport coefficients in underdoped cuprates. 

VI. GENERALIZED APPROACH TO TRANSPORT: 
T-MATRIX THEORY BELOW T c 

We turn now to transport properties below T c , beginning 
with the superfluid density. It should be clear from Figs. ^ 
and|2]that the order parameter and the excitation gap are to be 
distinguished in the superconducting state. We can thus say 
that the same T c as calculated in Section||I]via the self consis- 
tency conditions at the (extended) Hartree-GL level will show 
up, for example in p s {T). In particular the superfluid density 




FIG. 4: TDGL coefficients a' and <jq vs. impurity concentration 
m, for several values of scattering strength u from the unitary to the 
Born limit, in a self-consistent T-matrix treatment (Ref. 25), where 
the nonmagnetic, s-wave scattering potential is assumed, as given by 
u(x) = u6(x). 



must necessarily couple to the pair fluctuations in the super- 
conducting state in such a fashion that p s (T) reflects the su- 
perconducting order parameter, rather than the fermionic exci- 
tation gap. This coupling of pair fluctuations to p s (T) can be 
contrasted with the way in which collective (phase) mode con- 
tributions enter into p s (T) at the BCS level. These terms are 
only required to preserve general gauge invariance and these 
bosons do not affect p s when it is computed in the transverse 
gauge. The pair fluctuations we discuss here are necessary for 
a consistent calculation of the superfluid density, even in the 
transverse gauge. 

In this section we decompose the transport contributions 
into two types of excitations of the condensate: fermionic and 
bosonic. It is well known 7 that bosonic contributions to trans- 
port coefficients in the (less self-consistent) Gaussian-TDGL 
theory of pairing fluctuations are associated with Aslamazov- 
Larkin diagrams. The lowest order T-matrix theory introduces 
additional diagrams called the "Maki-Thompson" and "Den- 
sity of States" contributions. These latter two may be viewed 
as fermionic contributions. 

Given the self energy and form of self consistent T-matrix, 
Ward identities can be used to characterize the fermionic and 
bosonic contributions to transport, at the extended Hartree 
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level. Again, Aslamazov-Larkin (interpreted as associated 
with non-condensed bosons) and Maki-Thompson diagrams 
(interpreted as associated with fermionic excitations) are 
present, but both contain bare as well as dressed Green's func- 
tions. The fact that the combination GGo appears plays an 
important role, throughout our discussion. This pair suscepti- 
bility is related to the usual Gor'kov F function. 

Because of the form of the pair susceptibility and the fact 
that below T c , /j pQjr = 0, we show below that the fermionic 
terms combine to yield the usual BCS contributions to trans- 
port, which now depend on the full gap A, as distinct from the 
order parameter, A sc . That the fermionic contributions con- 
spire to be of the BCS form may seem natural at one level, but 
on another level, this is highly non-trivial since at very strong 
coupling A is essentially T-independent below T c , as can be 
seen from Fig. [3] Thus, the fermionic contributions to trans- 
port are nearly T independent, in striking contrast to what is 
found in the BCS or weak coupling limit. To demonstrate 
how the bosonic and fermionic contributions enter, we begin 
with a formulation of the superfiuid density in the presence of 
general self energy effects. 



A. Superfiuid density: General formalism in transverse gauge 

The electromagnetic response in our extended Hartree the- 
ory is constructed from the two additive contributions to the 
self energy shown in Figure [5^. The superconducting self en- 
ergy £ sc contains an anomalous reversal of the fermion line 
due to the creation of a condensed pair. While Y, pg apparently 
also has this reversed line, the charge flow is compensated by 
the fermion pair propagator, t pg , and is thus a conventional 
diagram. The quantity t pg , in turn, depends on one full and 
one bare Green's function, each of which, as shown below, 
leads to additional (Aslamazov-Larkin like) diagrams in the 
electromagnetic response. 

We define the superfiuid density, n s , in terms of the mag- 
netic London penetration depth Al as 



A 



-2 _ A*oe n s 



(44) 



where /xq is the magnetic permittivity. For convenience, we 
will set fi = e = c = 1. Note on a lattice, n/m should be 
replaced by 



d 2 e k 
dk^dkv 



Uk , (ji = x, y, z) , (45) 



where is the fermion density distribution in momentum 
space. For the quasi-2D square lattice, the in-plane mass ten- 
sor is diagonal, and (n/m) xx = (n/m) vy . 

We consider the in-plane penetration depth which is ex- 
pressed in terms of the local (static) electromagnetic response 
kernel K(0) in linear response theory™, 



\f = K xx {0) = (-) -P XX (Q) 

\Tfl / xx 



(46) 



where K is defined by 

MQ) = p^MQ) ~ (-) MQ) = -k»u{Q)A v {Q) , 

\77l/ fxv 

(47) 

and the current-current correlation function 

P„AQ)= dre tn - T 0V(cbr)jV(-q,0)> (48) 
Jo 

= -2 A^ M (K, K + Q)G(K + Q)X V {K + Q, K)G(K) . 



Here the bare vertex is given by 

X(K, K + Q) = V k e k+q/2 = - (k + § 



(49) 



where, for simplicity, we have used the dispersion for jellium 
in the last step, although the generalization to the lattice is 
straightforward. The electromagnetic vertex can be written 
in terms of the corrections coming from the two self-energy 
components as 



A EM = A + SA pg + 5A SC , 



(50) 



where SA pg is the pseudogap term. For illustrative purposes 
we specialize to the case ip = 1 (as occurs in s-wave pairing). 
The superconducting vertex contribution is given by 



SA SC (K + Q,K) 
= A 2 SC G (-K - Q)G (-K)X(K,K + Q) 



At the Q — » limit, it satisfies 



-SA SC (K,K) = 



dk 



(51) 



(52) 



This equation should be contrasted with the T > T c Ward 
identity 



5A pg (K,K) = 



dk 



(53) 



The difference in sign between Eqs. d52l and d53l is fun- 
damental and arises precisely from the anomalous nature of 
the E sc diagram. Indeed, Eq. (I53> is equivalent to the state- 
ment that above T c the paramagnetic and diamagnetic current 
contributions to the Q = response precisely cancel. This 
cancellation appears in the superfiuid density calculation pre- 
sented in the next section. By contrast Eq. J52i expresses 
the failure of this cancellation for the superconducting com- 
ponent, which is naturally associated with a Meissner effect. 

The vertex SA pg may be decomposed into Maki-Thompson 
(MT) and two types of Aslamazov-Larkin (ALi, AL^) dia- 
grams, whose contribution to the response is shown here in 
Fig- 05. We write 



SApg = SA M t 



5A\ L + 5A\ L {A) 



(54) 



Note that the AL2 diagram is specific to the extended Hartree 
scheme, in which the field couples to the full G appearing in 
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the T-matrix through a vertex A. It is important to distinguish 
the vertex A from the full EM vertex of Eq ( I50> . In particular, 
we write 



A = A + SA pg — SA S 



(55) 



where the sign change of the superconducting term (relative 
to A EM ) is a direct reflection of the sign in Eq. d52i . 

We now note that for ip = 1 there is a precise cancellation 
between the MT and AL\ pseudogap diagrams at Q = 0. 
This follows directly from the Ward identity 

Q ■ \(K, K + Q) = Gq 1 {K) - Go 1 (K + Q) , (56) 

which implies 

Q ■ [8k\ L (K,K + Q)+ SAmt (K,K + Q)]=0 (57) 

so that 8A\ L {K, K) = —5Amt(K,K) is obtained exactly 
from the Q — » limit. 

Similarly, it can be shown that 



A(K, K + Q) = G-\K) - G-\K + Q) 



(58) 



The above result can be used to infer a relation analogous to 
Eq. (1571 for the AL2 diagram, leading to 



SA pg (K,K) = -SA M t(K,K) 



(59) 



which expresses this pseudogap contribution to the vertex 
entirely in terms of the Maki-Thompson diagram shown in 
Fig. I3J). It is evident that 8Amt is simply the pseudogap 
counterpart of 8A SC , satisfying 



-SA M t(K, K) 



dk 



Therefore, one observes that for T <T C 

8E pg (K) 



SA pa (K,K) 



dk 



(60) 



(61) 



which establishes that Eq. ( 15 3> applies to the superconduct- 
ing phase as well. As expected, there is no direct Meissner 
contribution associated with the pseudogap self-energy. 

B. Superfluid density: T-matrix approximation 

The final expression for the superfluid density is obtained 
by rewriting Eq. ( 145 \ by integration by parts, 



(m) Q /3 2 ^ dk a 3kp G ^ 2 ^dk a 
= -2^G (K)—\— + — 



9e k dG(K) 
9S„ 



K 



dk a \dkf3 dkp dtp 



(62) 



Note here the surface term vanishes in all cases. 



By inserting Eqs. (I62> and J48l > in Eq. J46l > one can see that 
the pseudogap contribution to n s /m drops out by virtue of 
Eq. (16 II . The in-plane superfluid is isotropic and is given by 



m 



K 



dk x 



SA SC (K, K) x 



8Z SC {K) 



dk x 



(63) 



Equation d63i can be readily evaluated using the su- 
perconducting vertex and the superconducting self-energy 
T, SC (K) = —Al c Go(—K)ip^ associated with our GGo-based 
T-matrix approach. In addition, we introduce an approxima- 
tion in our evaluation of G via Eq. ( 136b . to find 



m 



A 2 

2 E% 



dk, 



1 



- 2/CgjQ 
2E k 



+ f'(E k ) 



4 dk T dk T 



(64) 



This result includes the general ip^ factor for the k- 
dependence of the gap and order parameter, whereas it has 
been neglected above in Eqs. ( I51> -( I61> . where it substantially 
complicates the analysis. 

Note that in the absence of a pseudogap (i.e., at weak cou- 
pling), A sc = A. Then Eq. (I64> is just the usual BCS formula. 
More generally, we can define a relationship 



(rh\ _ A 2 SC rn s 
V m ) A 2 V m 



BCS 



(65) 



where {n s /m) BCS is just (n s /m) with the overall prefactor 
A^ c replaced with A 2 in Eq. J64i . Obviously, in the pseudo- 
gap phase, (n s /m) BCS does not vanish atT c . That the results 
are so similar to their BCS counterparts is due to our Hartree 
treatment of pairing fluctuations. 

Finally, from Eqs. d37i and d65i we may write 




BCS 



(66) 



Here the first term represents the contribution to n s from the 
usual fermions, albeit with an unusual T-dependence of the 
gap (see Fig. [3}- The second term indicates that n s is addi- 
tionally depressed by bosonic pair excitations which ensure 
that n s vanishes prematurely at T c , rather than T*. 

Finally, from Eq. i53\ we can infer that the self energy ap- 
proximation of Eq. (I36> implies an approximation on SA pg , 
so that 

SAmt (K, K) « A 2 pg G (-K )G (-K )X(K, K) . (67) 

We build on this internal consistency argument in what fol- 
lows away from Q = 0. 



C. Q electrodynamics 

The real, fl ^ part of the in-plane optical conductivity 
can be expressed as 



(68) 
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FIG. 5: Self energy contributions (a) and response diagrams for the 
vertex correction corresponding to E P9 (b). Heavy lines are for 
dressed, while light lines are for bare Green's functions. Wavy lines 
indicate t vq . 



which is related to the superfluid density through the f-sum 
rule 

^ + - j°° dSl = (-) . (69) 

m it J \mJ xx 

Just as for the superfluid density, the optical conductivity 
involves the same set of Maki-Thompson and Aslamazov- 
Larkin diagrams. In this section we will regroup terms so as 
to identify explicit fermionic and bosonic contributions 

A = A -f- 6 Ay ermions + SA bosons . (70) 

Similarly, it follows that the optical conductivity contains two 
contributions 

o-(Q) = a fermions (n) + a bosons {n) , (71) 

where a f ermlons comes from the A + 8 A f ermions portion of 
the vertex. 

It is not unreasonable to take Eq. J67> a step further and 
apply it to general Q, so that (below T c ) 

6A MT (K,K + Q) 
» A 2 pg G (-K-Q)G {-K)X(K + Q,K) . (72) 

In effect what this approximation is saying is that for 
fermionic degrees of freedom the bosons enter primarily as 
an excitation gap contribution. This approximation is justi- 
fied by the same reasoning which leads to Eq. d36i . using the 
divergence of t pg (Q) at Q — 0. 

From Eq. J72I and Eq. J52I it follows that the pseudogap 
and superconducting condensate terms add in a natural way to 
introduce the full excitation gap A: 

8 'A f ermions 

(K, K + Q) 
= SA MT (K,K + Q) + 5A SC (K,K + Q) 
= A 2 G (-K -Q)\(K + Q,K)G (-K) . (73) 



This term, combined with the "density of states" contribution 
from A, indicates that the fermionic contributions to general 
transport coefficients are to be calculated within a BCS frame- 
work, but with the full gap A, which is to be distinguished 
from the order parameter. Though these contributions to trans- 
port are formally similar to BCS theory, their T-dependence 
may differ considerably due to the weak T-dependence of A 
in the strong pseudogap regime. 

To characterize the direct contribution from bosonic de- 
grees of freedom to transport, which are associated with AL 
diagrams, we must treat their full dynamics. We define 

SA bosons (K,K + Q) 
= 6A AL (K,K + Q)+5A 2 AL (K,K + Q) (74) 

and turn now to the bosonic contribution to conductivity. 



D. TDGL approach to bosonic transport 

The integrated conductivity of the bosons can be deduced 
from the fact that it accounts for the difference between the 
true superfluid density and the BCS prediction, leading to the 
constraint 

or 00 , A 2 /n \BCS 

- dna bosons (n,T) = -^(^-) (T). (75) 

TT JO A V 777 / 

The bosons make a maximum contribution at T c . At this tem- 
perature, bosons can account for as much as 90% of the spec- 
tral weight at T c in a strongly pseudogapped superconductor 
with T* /T c = 10, while their contribution vanishes in the 
weak coupling limit T* — » T c . The boson weight vanishes at 
T = at all couplings as a consequence of the condensation 
of all bosonic pairs. This pairing fluctuation bosonic contri- 
bution is not limited to a narrow region near T c , but extends 
well into the superconducting state. 

We define the bosonic response Pboson as the contribution 
to P given by 6 Arsons- These terms each involve a pair of 
T-matrices, and to leading order in frequency, Pboson may be 
written as 

Pboson (Q) 

= -2J2[^al(K,K + Q) + SA al (K,K + Q)] 

K 

xG{K + Q)\{K + Q, K)G(K) 

w Y,MP,P)t P9 ( p + Q)MP,P)t pg (P) . (76) 
p 

Here A t is the vertex for t pg approximated as 

A t (P,P)«2^p (77) 

where we have used the T-matrix expansion of Eq. i4ll . The 
quantity A t = zA t (P, P) where z = 1 in the normal state 5 
and is modified in the superconducting state. Reasonable es- 
timates of z(T) below T c may be obtained from Eq. ( 1691 . 
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If we presume Eq. H6\ holds for a range of low frequen- 
cies we may infer a simple expression for the (in plane) ac 
conductivity 



bosons 



m . E (5)7 *> 

xi(p. E + ») t ' E '- t n ' E + "» , ,78, 



where we now use A(p, SY) = —2 Im t pg (p, ft + iO) for the 
bosonic spectral function. More generally, at higher frequen- 
cies the internal fermionic structure via the Q dependence 
of the boson vertex must be included. This structure results 
from the individual coupling of radiation to each constituent 
fermion in the pair. However, it is reasonable to assume that 
the compositeness of the pairs will not be resolved by radi- 
ation of wavelengths larger than the pair size or frequencies 
below the pair-breaking energy A. For q = and £1 < A, 
then, we argue the bosonic vertex functions are well approx- 
imated by the velocity p/M (or a constant multiple thereof.) 
Indeed, for calculations of the ac conductivity, Varlamo\s2 8 .has 
argued that an analogous approximation (for Gaussian fluc- 
tuations) is valid near T c where the pole structure of the T- 
matrix causes its frequency dependence to dominate that of 
the Green's functions at the vertices. At the extended Hartree 
level this pole structure is present at and below T c due to 
the vanishing of /J, pa i r . Relative to Gaussian theory, the Q- 
dependence of the boson vertices will be further suppressed 
for fi < A through the appearance of the gapped fermion 
propagator G in all vertex subdiagrams. 

The result derived above for the bosonic contribution to the 
ac conductivity is essentially the same as would obtain from 
"true" bosons, except for the constant factor Oq 2 . To see this 
we review the ac conductivity in a system of free bosons of 
charge e* , mass M*, and chemical potential //*, in contact 
with a reservoir of fermion pairs. This calculation 26 represents 
a generalization of standard TDGL-like schemes, away from 
the Bose condensation temperature. The simplest physical 
picture which allows for an exactly solvable conductivity in 
the presence of quantum dissipation assumes that the bosons 
interact with a reservoir of localized 26 fermion pairs (treated 
as having ideal gas Bose-Einstein statistics). This gives rise to 
a boson self-energy (f2) without introducing vertex correc- 
tions to the electromagnetic response. Such a model yields 26 
an ac conductivity given by 



p 

b(E) - b(E + H) 




f irp 

—A(p,E)A(p,E + Cl) 

Z7T 



(79) 



Here b(E) is the Bose statistical function and A = — 2Imi? is 
the boson spectral function. The boson propagator is given by 
^(q,^)- 1 = Cl-q 2 /2M* + fi* -E B (0,). The boson vertex 
here is the velocity p/M* . The boson self-energy arises from 
scattering processes into and out of the thermal reservoir. 

To compare with the pairing theory, we note that while the 
boson self-energy in the above model arises from scattering 



into the reservoir, the self-energy of bosons in the pairing 
model (whose imaginary part we may regard as Im^ 1 ) arises 
from pair dissociation and recombination processes. We note, 
however, that fundamental differences between fermion pairs 
and true bosons remain in the analytic structure of the respec- 
tive T-matrix and boson propagator B. For true bosons, the 
real and imaginary parts of obey Kramers-Kronig rela- 
tions and vanish in the high energy limit. The propagator B 
then reduces to its "bare" form (SI ~ q 2 /2M* + fi*)^ 1 . The 
T-matrix has the structure i~ 1 = g _1 + \ where the pairing 
susceptibility \ satisfies the same causality constraints. The 
vanishing of \ m the high-energy limit leaves the asymptote 



t 



I'll 



g and the T-matrix loses all energy and momentum 
structure due to the dissociation of all pairs. However, this 
difference is not expected to be relevant for conductivity cal- 
culations done below the pair-breaking energy scale. 

We now consider the evaluation of Eq. ( I78t using the ex- 
panded T-matrix of Eq. 14 H . neglecting the £1 2 term. Dissipa- 
tion in a at nonzero frequencies requires cl'q ^ 0, which below 
T c requires impurity scattering, as discussed in Appendix Icl 
Here, we focus on the behavior of the conductivity as a func- 
tion of the ratio v = Oq/oq, which in turn enters T(q, f2) in 
the expanded T-matrix approximation (Eq. (IB2» as vVt. We 
note here that in the TDGL approach for Gaussian fluctuations 
near T c , this ratio is typically large (of order (Ep/T c )) due to 
the fact that a' is a measure of particle-hole asymmetry. As 
shown earlier in Fig.|4] v tends in the extended Hartree theory 
to increase rapidly from zero, becoming of order unity with 
the introduction of a modest concentration of impurities. 

TDGL calculations are typically done in the classical 
regime O < T, which allows the simplification [b(E) — 
b{E + fi)]/n -> T/[E(E + SI)] in Eq. (EH). The conduc- 
tivity below T c is then found to have the characteristic S7 -1 / 2 
dependence 28 for all values of v: 



1/2 



(80) 



To compute the conductivity in the superconducting state up 
to £1 ~ A, however, requires the inclusion of quantum sta- 
tistical factors when A ^> T c , in which case the v parameter 
affects the frequency-dependence of /j boson . While the low- 
frequency limiting behavior is £7 -1 / 2 independent of v, a falls 
faster than £7 -1 / 2 outside the classical regime before crossing 
over to Q 1 / 2 behavior at higher frequencies. This crossover 
may extend over a large frequency range, depending on the 
value of v, so that generally a is characterized by £l~ x l 2 at 
low frequencies and a long high frequency tail which cannot 
be integrated to infinity. Since the integrated bosonic weight 
is finite, it is reasonable to expect that this expression for the 
conductivity is cut off above the pair-breaking scale (at sev- 
eral times A(0)), where the TDGL formulation is known to 
break down. 



VII. CONCLUSIONS 

This paper is based on the observation that the BCS ground 
state wavefunction has a more general validity 1,29 . By increas- 
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ing the strength of the attractive interactions, superconductiv- 
ity in this state progressively takes on the character of Bose 
Einstein condensation. That the same T — wavefunction 
can apply to a system where pair formation and pair con- 
densation are associated with different energy scales provides 
support for the notion that there exists a mean field theory of 
a more general nature than that of strict BCS theory. In the 
more general case the various energies of the BCS picture are 
no longer degenerate: T* ^ T c and A ^ A sc . 

In this paper we have shown that this generalization of BCS 
physics is to be found in a treatment 4 - 6 of pair fluctuation ef- 
fects at the Hartree level. This establishes that the excitation 
gap A(T C ) is non-zero, or equivalently that there is a pseu- 
dogap, even in the weak coupling limit, as was observed ex- 
perimentally many years age 3 . Going beyond this previous 
work and into the broken symmetry phase, a self consistent 
analysis leads to the condition that the pair chemical potential 
Mpair = for T < T c , and that this is equivalent to the state- 
ment that the excitation gap A(T), assumes the BCS value 
at and below T c . It follows from the fact that the order pa- 
rameter is necessarily zero at T c , that the excitation gap and 
the superconducting order parameter are distinguishable be- 
low T c . This behavior is schematically illustrated in Figs. ^ 
and0 

Because it exists in conjunction with fermionic degrees of 
freedom, the ideal Bose character found here supports true 
superconductivity. This is demonstrated here by calculations 
of the superfluid density. More generally, transport proper- 
ties contain two types of contributions from both fermionic 
and bosonic excitations. At finite temperatures, the fermionic 
terms are well approximated by the usual BCS contributions 
to transport but with a highly-non -BCS-like, often tempera- 
ture independent gap, A(T). 

At this extended Hartree level, fluctuations in the order pa- 
rameter A sc and in the excitation gap A represent distinct 
channels for bosonic effects in the superconducting state. It 
is more convenient to decompose these into "condensed" and 
"non-condensed" bosons. The latter are the focus of the 
present paper, and we refer to these as "pair fluctuations" 
as distinct from fluctuations of the order parameter. Order 
parameter fluctuations were discussed in the context of our 
extended Hartree theory in earlier work 20 . Non-condensed 
bosons are present above and below T c , but absent (like 
fermionic excitations) at T = 0. These long lived, metastable 
states, even above T c (albeit, below T*) are associated with 
our extended Hartree treatment which introduces a pseudogap 
in the fermionic spectrum. This "gap" then inhibits dissocia- 
tion of the bosons into fermions. 

The dynamics of these non-condensed bosons is reasonably 
described by time-dependent Ginzburg-Landau theory. While 
the microscopic character of this TDGL changes abruptly 
from above to below T c , this approach offers a very powerful 
techniqu o 26 i 30 for addressing the pair fluctuation or bosonic 
contributions to transport. In this paper we have applied 
these results to calculations of the ac conductivity. One could 
equally well address the thermal conductivity and we can an- 
ticipate some of the ensuing conclusions. When both bosonic 
and fermionic excitations are present (i.e., away from T = 0), 



we expect that the Wiedemann-Franz law is violated. Because 
of the associated soft energy scale, bosonic contributions to 
the thermal conductivity are considerably smaller than their 
counterparts for the electrical conductivity. In this way the 
thermal conductivity k in the superconducting state should 
be well approximated by considering only BCS- like contri- 
butions, but with the anomalous temperature dependent exci- 
tation gap shown in Fig. [3] Moreover, because the bosonic 
contributions are unimportant in k, we expect to recover the 
well-known universal result 31 for this property. This is not 
necessarily true for the electrical conductivity. 

In this paper we have not made reference to specific phys- 
ical systems where our mean field theory may have some ap- 
plicability. In addition to short coherence length superconduc- 
tors, (among these the high temperature cuprates), the present 
picture may also be relevant 32 ' 33,34 to fermionic superfluidity 
in atomic trap experiments. 



Acknowledgments 

This work was supported by NSF-MRSEC Grant No.DMR- 
0213745 (JS, AI and KL), NSF grant No. DMR0094981 and 
JHU-TIPAC (QC). We thank Shina Tan, I. Ussishkin and A. 
Varlamov for useful conversations. 



APPENDIX A: T-MATRIX FOR SMALL A AT WEAK 
COUPLING 

If we replace the dressed Green's function in Eq. dlOl by 
the bare Green's function Go, we obtain 



<J 



1 



1 



l + 9Xo(Q) N(0) e + aiaj+CiV 



The self energy T,(K) is then given by 



(Al) 



E(K) = t(Q)Go(Q - A> 2 _ q/2 (A2) 
Q 

To calculate the pair susceptibility to the first order in self 
energy, we write the dressed Green's function in Eq. dlOt as 
G = Gq + GoSGo, which gives the first order correction 

Sx(Q) = E °o(Q - K)G (KMK)G (K) (A3) 

K 

Since t(Q) is sharply peaked around Q = we disregard the 
momentum dependence of Gq in Eq. (IA2> . and obtain 



E(P) « -G (-K)A 2 ^ 



where 



A = — 



d 3 q 
(2ir) 3 



i(q, LO n = 0) 



(A4) 



(A5) 



Putting this back into Eq. ( IA3I we find 

S X (0) = -fro/3 2 A 2 = b /T 2 J ^|i(q, = 0) (A6) 
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where 6 /T 2 = Ek G 2 (K)G 2 (-K)^ 2 k = N(0)^ (this 
last equality holds only for s-wave pairing). Using the expres- 
sion for the T-matrix corrected by #x(0) via Eq. H2\ 



t(Q) 



l 



N(0) e + a\n n \+gq*-6x(0)/N(0y 
the self-consistency equation reads 

b Q f d 3 q 1 



(A7) 



S X (0) = - 



N(0)Tj (2tt)3 e + gq2-5 X (0)/N(0) 



(A8) 



where 



2M 



(B5) 



is quadratic. This defines the effective pair mass, M . Below 
T c , we have 



fJ-pair(T) =0, T <T C 



(B6) 



Using Eq. JB U this condition can be shown to be equivalent 
to Eq. d38t . the familiar BCS gap equation. 



APPENDIX B: T-MATRIX FOR ARBITRARY A, BELOW T c 
FOR GENERAL COUPLING 



APPENDIX C: IMAGINARY PART OF THE INVERSE 
T-MATRIX - BOSONIC LIFETIME NEAR T c AND BELOW 



FromEq. (36) and G (P) = l/(iw„ -e p ) we find G{P) = 
(iuj n + e p ) /[(iuj n ) 2 — Ep]. Here e p is the electron normal 
state dispersion measured with respect to the chemical poten- 
tial (j,, while E p — y^ep + A 2 ip 2 . Thus after performing the 

Matsubara sum and analytical continuation ifi n — > SI + i0 + , 
Eq. (1 1 Oft becomes 



x(q,«) = E 



1 - /(£k) - /(Ck-q) 



e k-q 
/(£*) " /(£k-q 



- ^0^ 



e k _ q + SI + i(U 



Vk-q/2 • (Bl) 



In the long wavelength, low frequency limit, one can ex- 
pand the inverse T matrix as: 



C(q,^) 



9 M 
ctiSl 2 



x(<L«) 
■f a^Sl - 



2M 



Impair 



iT. 



q,Oj 



(B2) 



The linear contribution (q) is absent due to the inversion 
symmetry (q <-> — q) of the system. 

In the weak coupling limit, the ratio a' Q /a\ is vanishingly 
small; when the system has exact particle-hole symmetry 
(e.g., a 2D tight binding band at half filling with a nearest 
neighbor hopping), a' Q vanishes. In this case the dispersion 
determined via 



t" 1 (q,r!)=0 



(B3) 



is linear in q, Sl q ~ cq, which shows up in the dispersion, 
only in the very weak coupling limit or where there is exact 
particle-hole symmetry. 

In the absence of particle-hole symmetry, as g increases, 
a' /ai increases, thus a' fl gradually dominates and we find 
the important result: S7 q ~ q 2 . For any finite g and arbitrar- 
ily small q, the dispersion is always quadratic, at the lowest 
energies. 

We are interested in the moderate and strong coupling 
cases, where we can drop the aiSl 2 term in Eq. iB2\ . and 
hence we have 



si-sx, 



■iT. 



(B4) 



In this Appendix, for definiteness, we assume d— wave pair- 
ing so that = 2cos(20), where <\> is the polar angle. 
By taking the imaginary part of the inverse T-matrix using 
Eq. ( IB 1 1 we obtain the expression 

r q ,o = V -f( E ^ -/(ek-q)]«k*(3k + £k- q - n) 

CLf\ _ L 



lf(E k ) -/(ek- q )]wk*(^- £k- q + n) 



^k-q/2 • 

(CI) 



For small £1 <C T, and setting q = we can expand the Fermi 
functions to first order in SI: 

r ,n = -4^E \f'(E k )ul5(E k + e k -n) 

(In _ L 



■ f'(E k )vl6(E k - e k + fi) 



(C2) 



Here T(0, S7) reflects the rate of decay of non-condensed 
bosons into a bare and dressed fermion. The 5 functions in 
the above expression determine the energies of these fermions 



Ek and E k to be of order 



At temperatures of interest 



(T < T c and T > T c ), away from the nodes (where tp k = 0) 
and for SI < A this energy is large. The only appreciable con- 
tribution then comes from parts of the Fermi surface near the 

nodes (of angular width of the order of y^\), giving a higher 
power than linear in £1 dependence of T(0, SI) and <Zq = 0. 

This result is improved upon by including the effects of 
fermion-boson scattering. Quite generally in the absence of 
impurity scattering one finds that 

<' = -jE^(k,-e k )/'(e k ) (C3) 



where AQt, E) =-2 Im G(k, iuj n — > E + i0 + ) is the fermion 
spectral function. Fermion-boson scattering broadens the 
spectral function and A(k, —e k )f'(e k ) is generally non-zero 
in the normal state for all k, implying that cl'q ^ 0. One way 
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to implement this is to introduce an inverse lifetime 7 into the 
expression for the pseudogap self-energy: 



A 2 2 



LO + 6k UJ + £k + «7 



(C4) 



Below T c however A(k, — £k) in Eq. dC3> vanishes due to 
the sharpness of the superconducting self-energy, as is easily 
checked by inspecting the expression for the spectral function 
below T c following from the self-energy Eq. iC4\ 



leading to the following expression for the spectral function 
A(k, u), above T c 



(cc 2 -£ 2 ) 2 +7 2 (co-e k ) 2 



(C5) 



As is evident from Eq. dC5> , the quasiparticle peaks are 
broadened as a result of a non-zero 7. Thus, above T c the 
imaginary part of the inverse T-matrix (or, equivalently, of the 
pair susceptibility) is now 



Imx(q, O) = r(q, fi) = 7A 

':<). ~. 1 2 .. . P" 



2 

pg 



X 



„2 ,„2 



1 - /(» - gk-q) - /(fk-q) 

^ [(O- e k - q ) 2 - i^] 2 + 7 2 (n - e k - q - e k ) 2 ^^-4/2 



Since the quasiparticle peaks are broadened, the boson can 
now decay into states near the Fermi surface not only near the 
nodes but everywhere else on the Fermi surface. r(0, 51) is 
now linear in SI and the coefficient of proportionality is easily 
found to be 



d_ 

on 



lmt -1 (0,fi) 



7A 2 



pg \ " 



n=o 



T 



4 7 2 



(C7) 

The second term in the denominator of the summand in the 
above expression can be safely neglected if 7 is taken to be 
much smaller than A pg , since the main contribution to the 
sum comes from the vicinity of the Fermi surface (where e k 
is small). The only place where this term is not negligible in 
comparison with Ap g ip^ is in a small region near the nodes 

(angular width oc -7^-). However, this correction can also be 
neglected, and therefore 



// 

a 



ATA 2 

pg 



E 



cosh 



(C8) 



2A 2 v 1 7 (uj + e k ) 



(co + e k ) 2 (^ 2 -£ 2 ) 2 +7 2 (w 2 



(C6) 



I - A 2 C ^) ■ 

(CIO) 

Thus we can infer that <Zq = 0. Nevertheless, a small amount 
of impurity scattering will restore a non-vanishing g'q, since it 
is likely to produce a non-zero A(k, — e k ) in Eq. dC3i . In the 
d— wave case, the impurity renormalized imaginary part of the 
pair susceptibility has the following form^ 

/OO J 

x A (n - u , q - k)[f(Q -w)- f(u>)]<pl_ n/2 , (CI 1) 



where both lmG R (the spectral function of the full Green's 
function G) and A$ (the spectral function of the bare Green's 
function Go) are dressed with the impurity self-energy. To see 
the order of magnitude of this effect, we make a crude ap- 
proximation of this formula by assuming that only the full 
G is dressed by impurities but not Go which is equivalent 
to the statement that expression dC3t with renormalized A is 
valid. We model the impurity self-energy in the Born limit 
by Sj TOp (w) = — is\u>\ (as is reasonable close to the Fermi 
surface in the d— wave case) , where s is a dimensionless con- 
stant parameterizing the concentration of scatterers. Assum- 
ing s«l (clean limit), we obtain 



A(k,-e k ) = -2 



s|e k |A 2 ^|_ 

2 C 4 



AM" 



4s 2 e 



(C12) 



2T 



which, in conjunction with Eq. 1C3I and identifying — s|e k | 
with 7 yields an expression identical to Eq. (IC7b . Using rea- 
soning similar to that leading to Eq. dC9> , we arrive at the 
following estimate for a ' in the presence of impurities: 



which, assuming that T <C /i, after a standard procedure 
gives, as expected, a non-zero result 



a'o = 2 In 2 iV(0) 



sT 
A 2 



(C13) 



iV(0)7 

A 2 

pg 



(C9) More complete numerical results are shown in Fig.|4] 
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